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We review the difficulties of the generalized Chaplygin gas model to fit observational 
data, due to the tension between background and pcrturbative tests. We argue that such 
issues may be circumvented by means of a self-interacting scalar field representation 
of the model. However, this proposal seems to be successful only if the self-interacting 
scalar field has a non-canonical form. The latter can be implemented in Rastall's theory of 
gravity, which is based on a modification of the usual matter conservation law. We show 
that, besides its application to the generalized Chaplygin gas model, other cosmological 
models based on Rastall's theory have many interesting and unexpected new features. 

PACS numbers: 04.50.Kd, 95.35.+d, 95.36.+X, 98.80.-k 



Large scales observations indicate that about 95% of the matter content is com- 
posed of two exotic components, one (dark matter) with zero effective pressure and 
responsible for the formation of local structures, and the other (dark energy) with 
negative pressure and responsible for the present stage of accelerated expansion. 
They are exotic in the sense that no elementary particle of the standard model 
is able to represent them. Indeed, they do not emit any kind of electromagnetic 
radiation and up to now were only indirectly detected through their gravitational 
effects. 

The conventional cosmological model contains the usual known components 
(baryons, radiation and neutrinos, understood as all the particles forming the stan- 
dard model), supplemented by dark matter and dark energy. The former is rep- 
resented by a pressureless component that decouple from radiation much earlier 
than the baryonsf^ whereas dark energy is modelled by a cosmological constant A, 
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allegedly representing vacuum energy associated to quantum fields!^ The standard 
model of cosmology is named ACDM. 

The most important sign of the success of the ACDM is its ability to nicely 
fit both kinematic (i.e. regarding the background evolution history) and dynami- 
cal (i.e. regarding perturbations about a homogeneous and isotropic background) 
cosmological observations. The free parameters are essentially three: Hq, i.e. the 
value of the Hubble parameter today, f2 m o and f^Ao, i.e. the dark matter and dark 
energy densities parameters. Under the condition that the universe is spatially flat, 
i.e f2 m o + S^ao = 1) the number of free parameters is reduced to two. Recent analysis 
of CMB observations, 3 obtain H = 70.41^ km s" 1 MpcT 1 , n A0 = 0.728±g;gi| 
and fi m Q = 0.227 ± 0.014. Moreover, the baryon's fraction to the critical density is 
given by ftbo = 0.0456 ± 0.0016 and the curvature parameter of the spatial section 
is —0.0133 < Oko < 0.0084, indicating an almost flat universe. 

In spite of its success, the ACDM model faces many difficulties, both theoretical 
and observational. The first one concerns the nature of the components of the dark 
sector. The most quoted candidates to represent dark matter are neutralinos (stable 
supersymmetric particles) and axions (predicted by the Grand Unified models). 
No one of the underlying fundamental theory for these candidates has received 
experimental evidence until now, thus they remain theoretical proposals. In some 
sense, evidences for the existence of dark matter may be considered as an indirect 
proof that such fundamental theories may have something to do with nature. On 
the other hand, the dark energy component, the cosmological constant, is usually 
interpreted as vacuum energy emerging from quantum fields in curved spacetime. 
However, theoretical estimates lead to values up to 120 orders of magnitude larger 
than its observed value. On the other hand, if considered as a pure geometric term, 
the observed value of the cosmological constant requires an extremely precise tuning 
of its value. 

Moreover, there are at least two other observational aspects that bring troubles 
to the ACDM model. The first one is the so-called cosmic coincidence and concerns 
the fact that the present fractional densities for dark matter and dark energy are of 
the same order. This is very unlikely to happen since dark energy and dark matter 
scale very differently with time. Second, predictions of the ACDM model for non- 
linear structures in the universe do not fit numerical simulations, mainly in what 
concerns the existence of an excess of power in the mass agglomeration at galactic 
scales 3 

In view of these problems, many alternatives to the ACDM model have been 
proposed!^ For example, quintessence and if -essence. However, both have their spe- 
cific problems, like the extreme small value of the scalar particle in the quintessence 
model (i.e. ~ 10 -33 eV) and the stability issues for the if-essence model. An 
interesting proposal is the interacting models for the dark sector, which consider 
the possibility of the decay of one component into the other, and may shed some 
light on the cosmic coincidence problem. On the other hand, they suffer from the 
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necessity of introducing phenomenological interaction terms, due to the absence of 
a fundamental model for the decaying process. At the same time, thermodynamical 
aspects concerning such interac tion processes lead to constraints that seem to be in 
contradiction with observations 

Among the alternatives evoked to describe the dark sector of the universe, one 
interesting proposal are the so-called unified models. Their prototype is the Chap- 
lygin gas (CG)p' which is inspired by string theory: it is a fluid able to mimic dark 
matter in the past and dark energy in the future, thus properly reproducing the 
expansion history of the universe. However, the CG suffers from some difficulties 
when it is tested against observation. This led t o a more general formulation, the 
so-called generalized Chaplygin gas (GCG)P^^ which possesses a new parameter 
a (a = 1 recovers the original CG). Unfortunately, a tension appears between the 
predictions of the model and variou s observational tests: the background ones (e.g. 
SNIa, BAO, H) favour a < q JHIIZI whe reas the per turbative ones require a > in 
order to have a positive speed of sound jEHHIHlS] 

In order to cope with this problem a possibility is to abandon the fluid descrip- 
tion for the GCG, using instead a self-interacting scalar field. Unfortunately, this 
cannot be done by using a canonical scalar field. We show in what follows that an 
interesting way out is to use a non-canonical self-interacting scalar field as suggested 
by Rastall's theory of gravity, an alternative to General Relativity where the impo- 
sition of a strict conservation of the energy-momentum tensor is relaxed!^ Besides 
the alternative description for the GCG model, Rastall's theory may lead to many 
interesting cosmological scenario, as we describe below. 

Consider the GCG equation of state: 

where p is the density and a and A are positive parameters. Inserting this expres- 
sion in the continuity equation in a Friedmann-Lemaitre- Robertson- Walker (FLRW) 
metric, i.e. 

P + 3% + p)=0, (2) 

where a is the scale factor and the dot indicates derivative with respect to the 
cosmic time, one finds the following expression: 

p = p [A+(l-A)a-^+ a '>] Th , (3) 

where A = A/ 'pQ a and po is the present-time GCG density (the present time is 
fixed by requiring ag = 1). 

The tension between background and perturbative constraints can be easily 
noticed from the equation of state ([lj: while in principle such definition allows for 
any real value of a, at the perturbative level we have to deal with the square of the 
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speed of sound, which has the following form: 



„2 



dp aA 



dp p 1+a 



(4) 



which is positive only for a > 0. 

The source of t he t ension comes from the fact that background tests favour 
negative values of ajJJJ which are excluded at perturbative level due to the occur- 
rence of an imaginary speed of sound. This fact is due to the fluid representation 
adopted. It must be noted that such fluid representation emerges naturally from the 
Nambu-Goto action in the original CG formulation!^ The DBI scalar field formu- 
lation of the Chaplygin gas is completely equivalent to the fluid formulation, even 
at the perturbative level] ' ' On the other hand, the GCG has not the same clear 
connection with string theory (even if a DBI-type action can be written presenting 
similar properties to the strict DBI case). In consequence, there is no fundamental 
reason to stay in the context of string theory. Therefore, it is possible to conceive a 
more fundamental representation for the GCG. A first possibility is to interpret the 
CGC as a self-interacting scalar field. However, a canonical self-interacting scalar 
field possesses a speed of sound equal to the speed of light, as we prove below, 
limiting the possibility to use such description for a unified model of dark matter 
and dark energy. 

The adiabatic speed of sound of a generic if-essence scalar field has the following 
expressionffSl 

cl = **■ (5) 
Px 

where x is the kinetic term. A canonical self-interacting scalar field has the following 
expressions for the pressure and the density in the FLRW background: 

p* = Y +v W> P* = y-^M. ( 6 ) 

and it obeys the Klein-Gordon equation, which, in the FLRW background, reads 

4> + 3H^=-V^) , (7) 

where V^(4>) = dV/d<fi. Hence, the adiabatic speed of sound for a canonical self- 
interacting scalar field is c s = 1. This proves that a canonical scalar field cannot 
play the role of dark matter. 

So, how a self-interacting scalar field could represent the dark matter compo- 
nent? We have to consider a non-canonical scalar field. We adopt Rastall's gravity 
theoryjl^ where the energy-momentum tensor T^ v does not obey the usual conser- 
vation law. Instead, the divergence of T^ v reads 



(8) 
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where R is the Ricci scalar and k is a constant. The field equations in Rastall's 
theory read 

R^u - ^g^R = 8ttG (t^ - ^-^-g^T^j , (9) 

3*% = r ^—^—T'' l/ , (10) 

where 7 is a dimensionless constant connected to k. When 7 = 1, general relativity 
is recovered. Using the canonical form for the energy-momentum tensor of a scalar 
field, i.e. 

Tp„ = (f),^<t>, v - ]^g^(f) tP (f)- p + g^V{<p) , (11) 
we obtain the following coupled equations: 

R pv - \g„uR = <i>,„<t>,u - ^pW,«0' Q + g^d - 2 7 )u(</>) , (12) 

+ (3 - 2 7 )^ = (1 - i f p fjr ■ (13) 

From Eq. (|12l) , the following effective energy-momentum tensor can be read off: 



implying the following expressions for the energy density and pressure: 



pf = |^ + (3- 2 1 )V{<f) , pf = ^-^, 2 -(3-2 7 )U(0). (15) 

Using this expression in Eq. ((5} to evaluate the speed of sound, one finds 

2 _ 2-7 



7 



(16) 



This implies a vanishing speed of sound for 7 = 2. In this case, the non-canonical 
self-interacting scalar field based on Rastall's theory may represent dark matter. 
On the other side, from the non-perturbative point of view, it can represent dark 
energy by a suitable choice of the potential. We will explore this possibility in what 
follows. 

Let us consider a self-interacting scalar field, with the effective energy- 
momentum tensor (1141) . with 7 = 2: 

= <t>A,» - g^v{<j>) . (17) 

Inserting the FLRW metric, we have the following density and pressure associated 
with this scalar field: 

Pc , = <p 2 -V(<f>) , p$ = -V(<f>). (18) 
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Let us suppose that this density and pressure reproduce the background behaviour 
of the GCG model. Hence, in this case, we have, 



V(a) = 3n c0 Ag(a)- a/{1+a) , 



(19) 
(20) 



where g(a) = A + (1 — A) a 3 ( 1+Q ). Hence, in order to have a zero speed of sound, 
the scalar model must obey the following equations: 

1 



R 



= 0, 



(21) 
(22) 



where, just for future convenience, we have made the redefinition V(<fi) — » —V((f>). 

Let us inspect now the perturbative behaviour of this system, computing scalar 
perturbations in the density contrast. The perturbed equations in the synchronous 



coordinate condition read 



20 



S + 2*S-^6 = j>*-V4* 
a 2 a J 

d • fk 2 \ 
2* + 3-* + — + V H * = 4> 



(23) 
(24) 



where = 6<f> and 5 is the density contrast of the matter component. Using the 
scale factor as independent variable, the above system of equations takes on the 
following form: 



6" 

2$" 



a 



f'(o) 
f'(a) 



i 3 
6- - 



O 



2a3/ 2 (a) 



8 = <p 



a }{a 



a 2 P{a) + P(a) 



V 6 



* = cp'5' , 



(25) 
(26) 



where /(a) = a = ^/f2 m oa _1 + £l c {a)a 2 and O c (a) = fl c og(a) 1 ^ 1+a K 

Using a Bayesian analysis and comparing the theoretical predictions of our model 
with the 2dFRGS data for the power spectrum of matter distribution in the universe, 
we find a significant probability region for a < 0. The results are shown in figure 1, 
for the unification scenario where the dark sector is represented by the GCG only.^ 
Hence, if the GCG model is represented by a non-canonical scalar field, like the 
one suggested by Rastall's theory of gravity, the observational tension that plagues 
the GCG fluid model may disappear or, at least, be considerably alleviated. This fact 
may open new perspectives for the dark matter-dark energy unification program. 

However, this non-canonical scalar formulation obeying the Rastall's conserva- 
tion equation has some striking peculiarities. Let us consider, for example, scalar 
perturbations in the Newtonian gauge: 



ds 2 = a 2 (l - 2$)d?/ 2 - a 2 (l + 2<f>)~/ ij dx i dx i , 
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2 4 6 



Fig. 1. PDF for the GCG parameter a for the unification scenario, i.e., SldmO — 0- 



(27) 
(28) 

(29) 



The perturbed RastalPs equations (fT2 ^ -(fT3 |) becomes: 

V 2 $ - 3H {U<S> + $') + 7 (K 2 - W) * 
AttG [70q<50' + (3 - 2 1 )a 2 V^8(f\ 
U<S>.i + & t = AnG<t>' 8<t>, 
+ + (2W + H 2 ) $ + (2 - 7) (H 2 - W) $ 
4ttG [(2 - -f)4>' 5(j)' - (3 - 2~f)a 2 V,<f,5<f 

These equations can be combined in different ways in order to obtain a single 
equation for the potential $. These different combinations lead to the following 
equations: 

$" + 3W + (2W + H T 

(2n' + n 2 ) 

.2 



7 



7 



[-fc 2 $-3"H (%$ + $')] 



(3 



7^ 
$" + 3H$' 



-//■-<!» - :!7/^^ + $') - ^^(3 - 2 7 ) (H$ 



)$ = 
$') , 
$ = 
$') . 



(30) 



(31) 



7$> 

These equations are identical only in two cases: either 7 = 1 or if the perturbed 
quantities are not space dependent, leading to k — 0. In the first case, the Rastall's 
theory reduces to general relativity. In the second case we have a pure redefinition 
of the background quantities, remaining with a homogeneous universe!^ 

A way to give sense to perturbations in the scalar formulation of the Rastall's 
theory without reducing it to general relativity, is to consider, besides the scalar 
field, a hydrodynamic fluid. In such situation, the problem presented above disap- 
pear and perturbations can be consistently defined for any value of 7. Hence, this 
non-canonical scalar field requires matter in order to make sense perturbatively. 



March 3, 2013 16:25 WSPC/INSTRUCTION FILE textoV2 



8 Fabris et al 

This is not the only surprise that results from the new conservation law (|10[) . In 
some situations, the Rastall's cosmology can lead to the same achievements of the 
ACDM model, but with new features at the non-linear regime where the ACDM 
model faces some difficulties!^ 

Let us consider a two fluid model, with one of the components obeying the 
traditional conservation law. In this case, equation (jlOl) splits into two equations: 

Tr Ui = > r r ;M = o , (32) 

where the subscripts m and x indicate the matter and dark energy components. 
Choosing p m = et p x = —p x , we find the following equations, considering also 
equation 



H 2 



(3 - 2j)p x + 7 2 + 3 p„ 



(33) 



7-1 

p m + 3Hp m = , (3 - 2j)p x = —-—p m ■ (34) 
Integrating for the densities, we obtain, 

_ PmO_ _ PxO 1 - 1 ,„_>. 

Pm - fl3 , Px - 3 _ ^ + 2(3 _ ^Pm • {^) 

Inserting these expressions in the modified Friedmann equation, it results 

H = — (PxO+Pm), (36) 

that is, Friedmann equation for the ACDM model. In spite of this, the dark energy 
component is now dynamical and scales with time as 

_ PxO 7 - 1 /o-x 

P ^ 3 - 2 7 2(3 - 2 7 ) Pm ' [6{) 

This means that all the successes of the ACDM model in fitting the so-called kine- 
matic observation tests (Supernova type la, BAO, H(z)), are recovered by such 
Rastall's cosmology. 

Astonishingly, the ACDM structure with a dynamical dark energy component is 
preserved at perturbative level. In fact, considering linear perturbations (for simplic- 
ity, we consider now the synchronous gauge condition), after standard calculations, 
we end up with a single equation for the matter density contrast S m = Sp m /p m : 

S m + 2-S m - Ai:Gp m 5 m = . (38) 
a 

It is again the same equation of the ACDM model. Hence, in what concerns the dy- 
namical cosmological tests, requiring perturbative analysis (like the CMB anisotropy 
and matter power spectrum) , Rastall's theory is perfectly equivalent to the ACDM 
model. 
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But, there exist fluctuations in the dark energy component given by, 

5px = 2(3- 2 7 ) 5pm ■ {39) 

In order words, now dark energy agglomerates. Such structure is preserved e ven at 
second order perturbations, when neglecting vector and tensor contributions P2l 

This agglomeration of the dark energy component should have important con- 
sequences at non-linear level. For example, let us consider the spherical collapse 
model. The equation describing the evolution of a spherical perturbation is given 

by, 



Hi 



n p (u)- + i - n p (u) , (40) 



where is the scale factor computed at some initial time U at which the collapse 
begins and fl p = 1 + S is the density parameter of the collapsing region, determined 
by the density contrast of the fluid there contained. In the Rastall's cosmology 
described above, we have 

n p = l+ 2 -^p±6 m + l-^5 x , (41) 
5 — 37 5 — 37 

where the perturbations on the dark energy component appear. This term is, of 
course, absent in the corresponding expression for the ACDM model, as can be seen 
by setting 7 = 1 in the relation above. 

All the features described here involving the Rastall's cosmology, reveal the 
richness of such a proposal based on the modification of the usual conservation law. 
However, it is clear that a deeper analysis is required, possibly in connection with 
quantum effects in the universe. At the same time, the implications of Rastall's cos- 
mology for the formation and properties of non-linear structures is a very promising 
research program. 
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